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SUMMARY 


An analysis is presented for the deformation of a doubly curved 
thin plate under edge loads or surface loads for small deflections. 

This problem is approached from thin-shell theory so that the plate is 
to form part of a shell of revolution. The method developed is particu- 
larly useful for a plate whose radius of curvature in one direction is 
large compared with its length and width dimensions. The solution con- 
sists of an expansion about a parameter which depends on this fact. 

An analytical solution is presented completely for a plate with an 
arbitrary meridian curve of small curvature and loaded by normal edge 
loads on one pair of opposite edges. Numerical calculations for the 
deflection and moment distribution are presented for a particular 
meridian curve. For the meridian curve chosen for the numerical example, 
part of the surface had a negative Gaussian curvature. Results show that 
the deflections and bending moment are largest at the part of the plate 
with negative Gaussian curvature. 

The method is developed to the point that it may be applied readily 
to other problems of the deformations of doubly curved thin plates under 
edge or surface loads. The theory, however, is limited to small deflec- 
tions of the plate or shell considered. 


INTRODUCTION 


This report is concerned with the behavior of a doubly curved thin 
plate under edge loads or surface loads. This problem is considered in 
the following way. The doubly curved plate is to form a part of a shell 
of revolution bounded by two meridians and two parallels. The meridian 
curve is assumed to have a radius of curvature much larger than the 
radius of curvature of a parallel. This situation is clearly presented in 
most airplane fuselage panels. The introduction of a small parameter 
dependent on this fact allows the equations for equilibrium of the shell 
of revolution to reduce to ones with constant coefficients. The solution 
of this sequence of problems then leads to a complete solution of the 
problem. It may be noted that the method so developed is equally valid, 
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for analysis of the deformation of a shell of revolution with a meridian 
curve of small curvature and loaded in any manner whatsoever. 

The thin-shell equations valid for small deformation of an arbitrary 
shell of revolution are first presented. These, in general, lead to a 
set of differential equations with variable.^oefficients. The pertur- 
bation method is then applied to reduce these equations to ones with 
constant coefficients. An analytic solution is presented completely for 
a plate with arbitrary meridian of small curvature loaded by normal edge 
loads on one pair of opposite edges and stress free on the other pair. 

In addition, all edges are free from shearing stresses. For this type of 
condition two infinite families of algebraic equations are obtained which 
are solved by retaining a finite number. Finally, numerical calculations 
for the deflection and moment distribution are given for a particular 
meridian curve. 

The author is indebted to Dr. D. J. Peery, Head of the Department 
of Aeronautical Engineering, for his work in the formulation of the 
proposal for this project and to Dr. J. A. Sauer, head of the Department 
of Engineering Mechanics, for his understanding cooperation. The author 
is also grateful to Miss A. Koo who checked the major portion of the ana- 
lytical work and to Mr. T. Khammash who performed the numerical 
calculations. 

This work was conducted in the Department of Aeronautical Engineering 
of The Pennsylvania State College under the sponsorship and with the 
financial assistance of the National Advisory Committee for Aeronautics. 


SYMBOLS 


A = Gh 3 /l2 
a 

B = Gh 

C = Eh/(l - v 2 ) 

D = Eh 3 / 12(l - v 2 ) 
E 

f(°), F^ 1 ), . . . 
f(z) 


constant defined by equation (10) 


Young's modulus 

stress functions 

equation of meridian curve 
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* 

G 


shear modulus 

* 

G, H, K 


known functions given from solution of each 
approximation 


h 


thickness of shell 


l, 6 0 


define length and width of plate 


M g> Mg, M|g 


moment resultants defined by equations (3) 


N 0> W |0 


force resultants defined by equations (3) 


P £ 


surface l’oad in direction normal to middle surface 


R 0 , R| 


principal radii of curvature of middle surface of 
shell 


u, v, w 


displacements in 0, and {; directions, 

respectively 

ft 

x, y, z 


rectangular Cartesian coordinates 


2 2 
or, r 


fundamental magnitudes of first order of middle 
surface 


a n = nrt/2 




Pm = “V 2 




7 m = mrtb/2 




5 n = nrtl/2b 


; 


e 0> e |> ^0? 


extensional strains and shear strain given by 
equations (8) 


5 


coordinate normal to middle surface 




curvilinear coordinates on middle surface of shell 


V V T 


change-of-curvature functions given by equations (8) 


P 


small parameter defined by equation (10) 


P m = mn / 2 i 
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If 


V 

V n = nrt/2b 

T le> T ei 

T |£> T ee 

t 

2 


Poisson’s ratio 


shear stresses 
normal stresses 

function defined, by equations (22) 
"equivalent" potential function 


BENDING OF SHELLS OF REVOLUTION 


The differential equations governing the small deflection of shells 
of revolution are considered herein. In this section the equations 
pertinent to the later investigations are outlined. 


Coordinate System on Shell . 

Consider a shell of revolution located as shown in figure 1; in 
addition to the x, y, and z Cartesian coordinate system a set of 
orthogonal curvilinear coordinates 0, | , and ■ £ is chOBen on the 

middle surface of the shell such that the 0 and £ lines are lines 
of curvature of the middle surface and £ is the distance normal to 
the middle surface. In addition, the 0 and | lines are the parallels 
and meridians of the middle surface (for a shell of revolution). Note 
that £ is a parameter along the meridian plane and 0 is the polar 
angle measured from the xz-plane. The element of arc length will be: 

ds 2 = a 2 ^l + d| 2 + r 2 ^l + JjL^ d0 2 + d£ 2 (1) 

where a 2 = + ^2^ and Eg and R Q are the principal radii of 

curvature of the middle surface. It is noted that instead of | one can 
use z as the independent variable. Then the element of arc length in 
the | direction becomes 

ad5 -\/ 1 + (Hf dz (2) 


One also notes that r = r(z) is the equation of the meridian curve. 
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The sign convention as given, for example, "by Sokolnikoff and 
Specht (reference 1) is used. Thus if Tj_j is the stress tensor, then 

the subscript i indicates direction of normal to plane under consider- 
ation and j, the direction of the component of stress. Tensile stresses 
and compressive stresses are positive and negative, respectively. Shear 
stresses on a particular plane are positive if the normal forces are 
positive on that plane and if they are acting in the direction of posi- 
tive coordinate axes j otherwise, they are negative. 

The stress resultants and moments, as defined later, have the fol- 
lowing sign convention (see fig. 2). The resultant bending moments are 
positive if they cause positive stress on the positive side (£) of the 
middle surface. Resultant forces have the same sign convention as the 
stresses. 


Assumptions in Analysis 

The assumptions inherent in this analysis are the usual ones for 
thin-shell theory. They are amply presented in reference 2. Summa- 
rizing, they are: 

(1) Material is isotropic and follows Hooke's law 

(2) Thickness of shell is small compared with smallest radius of 

curvature of middle surface 

(3) Displacements are small compared with thickness 

(J+) Straight lines normal to undeformed middle surface remain 
straight and normal to deformed middle surface 

In addition, two other assumptions are made here. These are: 

First, the effect of transverse shear in the resultant force equations 
is small and can be neglected and, secondly, the effect of the displace- 
ments tangential to the middle surface of the shell in the changes of 
curvature is of higher order than that of the displacement normal to 
the surface. These assumptions are discussed more fully later in the 
text. 


Differential Equations of Equilibrium 

Force and moment resultants .- The stress resultants and moments are 
defined on the coordinate curves on the middle surface of the shell. 

They are given in units of force or moment per unit length of the middle 
surface. Thus the following definitions are made: 




with the integrals evaluated, between the limits -h/2 to h/2. 

The differential equations of equilibrium for the force and moment 
resultants may be obtained from physical considerations (e.g . , refer- 
ence 2) or by integration across the thickness of the shell of the 
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differential equations for the stresses (reference 3). Thus the 
differential equations for the stress resultants are: 


where 


&H) + &(*■*>) + If N « + S 


81 


84 


= 0 


CM 


+ - If + if % - 0 


(4b) 


-^ H + ar , r ^ e )--0 


(4c) 


p 4 = 


1 + 4)! 1 + 


h/2 
- h/2 


In these equations it has been stipulated that there are no tangential 
forces applied to the middle surface of the shell and no body forces 
are presented. The equations for moment equilibrium are 


H^e) 

80 

dfos) 

84 


+ ^|e) - raQe + ll Meg = o 


+ §^06 ) - ^4 - If % = 0 


%4 - N 6e 



%4 

R 0 


(5a) 


(5b) 


(5c) 
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Equation (5c) is not obtained from integration of the stress- 
equilibrium equation in the £ direction. It is implied by the rela- 
tion of equality of cross shears 

T el = T ie , . 

that is, if equation (6) is rewritten in the form 

^ + - T ie^ * 4) “4 ^ + k) + t STt9 ( 1 + %) 

and integrated across the thickness of the shell, equation (5c) results. 

Relations between force and moment resultants and strains of middle 
surface . - The relations between the stress resultants and the strains of 
the middle surface are obtained from stress-strain relations and defi- 
nitions (3) utilizing the assumption that ^1 + « 1 and ^1 + ~ 1 

in equations (3). (See references 3, 4, and 5*) These relations are: 

Bg = C(6| + V 5 g) l 


N e _ c ( € e + Ve ‘ 


N Se " N e! _ B7 ie 


Mo = D(Kq + VKt 


M| = D(lt| + VK 0 ) 


= M 0 1 = At 


where A = B = hG, C = - and D = — -r. 

12 1 - V 2 ' 12(l - V 2) 
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The extensional strains and changes of curvature may he given in 
terms of the displacement of an arbitrary point in the middle surface of 
the shell. Let u, v, and v be the components of the displacement of 
a point on the middle surface where u is along the tangent to a parallel 
circle, v, along the tangent to the meridian curve, and w, along the 
normal to the surface. Thus the extensional strains and changes of curva- 
ture of the middle surface are 


e 0 


_ 1 8 u + _ 1 _ 8 r_ y + w_ 
r 80 ra R e 


n 


1 dv w 

a 81 B| 


y = u 8 r ... 1 8 u , 1 8 v 
7 0 | "ar a r £0 


= l/l 8u i 8 w \ 
B - 30 r ig 2) 


+ JL 8r /v 
or d| \Rg 


1 8 w \ 
a 81/ 


y (8) 


K. = — d /y _ 1 8 w\ 
I a S|\R| a 81/ 


T 


1 8 r (u 
ar 



1 8_ 

a 



In the latter part of this report the above set of equations will 
be applied to doubly curved plates which have been "cut" from a nearly 
circular cylindrical shell. Thus the principal radius of curvature Rg 

will be quite large. Moreover, the assumption (h/R) « 1 (r is the 
smaller of Rg and Rq) implies that in these problems u and v are 
second order small in Kg, Kg, and t. Finally, if terms of higher 

order than h^ in the values of Kg, Kg, and t are rejected one 
obtains in place of the last three relations in equations ( 8 ) the 
following equations: 
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k = 8^w + 1 8r 8w 

^ 80^ a^r 81 81 

Kt = _i _8/l 8w\ 

1 a 8*1*86/ 

t = ji 8r 8w 1 8 (l 8w \ l _8_/i 8w \ 
ar 2 8| 80 ~ ^ d£V r 80 J “ r a0V a 81/ 


(9) 


One notes that in equation (5c) 
the order of h^/ R | and h-^Rg 


the terms M^g^Ri and Mgg/Rg are of 
and can be rejected, thus leading to 


N |0 = N fi | 


agreeing with equations (7). A further simplification can be made by- 
rejecting the transverse shears Q@ and Q| in the force equations 


since they are of the order of h^/fig and h 3 /R|, respectively, while 
the remainder of the terms are at most of the order of h/R| and h/R g. 


ANALYSIS OF CURVED PLATES 


In the following sections the equations will be applied to the 
bending of doubly curved plates which form a part of a shell of revo- 
lution. The equations as given above are, in general, intractable to 
the problem at hand. However, in the application, for example, to 
curved panels of an airplane fuselage the radius of curvature of the 
meridian curve is many times the length and width of the panel. Hence a 
solution to the problem can be obtained by an expansion in terms of a 
parameter which will be small when the radius of curvature of the meridian 
curve is large. Then the solution will appear, for example, in w as 

w = w^°) + pw^) + + . . . 


where p is the parameter concerned. 
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To this purpose the equation of the meridian curve is represented 
by the relation 


r = a + pf(z) 


( 10 ) 


where a is a constant, p is the small parameter dependent on the 
maximum deviation of the meridian curve from the straight line r = a, 
and f(z) is the equation of the meridian curve. 

It is’ more convenient in the following sections to use the coordi- 
nate z instead of |. Thus one replaces I by z in all equations 
above. 


The radii of curvature R^ and R q are given by the relations 


_1_ 

R e 




r|l + (dr/dz) 2 ^ ^ 

(d^r/dz 2 ) 

[l + (dr/dz)^ 2 


> 


(ID 


In accordance with the first paragraph of this section all dependent 
variables will be expanded in power series of the parameter p. Then 
solutions of the different approximations are obtained. For snail values 
of the parameter p, probably two approximations are sufficient to lead 
to a complete solution of the problem. 

The expansions required are: 


> ( 12 ) 


= u(°) + pu^ 1 ) + p 2 ^ 2 ) + . . . 


(o) ( 1 ) 2 (2) 

= *e + n * e + n e e + • 


K e = + + f-* 2 ^ 2 ^ + • 




u 


(13) 
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% 

Mg 


(o) + pNg^ + vN[q^ + . . . 

(1*) 

^ + liMg^ + p 2 Mg ^ 2 ^ + . . . 

(15) 


and similar expansions are obtained for e z , €0Z , k z , t, N z , Ng z , 
M z , and M 0Z . 

The expansions of the various terms have been summarized in appen- 
dixes A and B. The superscripts correspond to the order of the 
approximation. 


Equations for First Approximation 

The first approximation yields essentially the same equations as 
that of a circular cylinder, as would be expected. One can obtain two 
simultaneous equations for a stress function F'°' and the normal displace 
ment w(°) from the pertinent equations in appendix B as follows: Let 


% (0) - , 
9 Sz 2 

K (o) 1 

z " a 2 ^ 2 : 

N (o) _ dV 0 ) ■ 

0Z " a 5© dz 


(16) 


where the superscript corresponds to the first approximation. Then it 
is seen that equations (B13) for force equilibrium are automatically 
satisfied. Moreover, a compatibility equation can be constructed from 
the strains (equation (B7)) and is 


a 2 , ( °) 

O eg 


_1_ 

.2 


5© 2 


h 2 y 

1 d ^ Q z 
a 8z b6 


1 8 2 w ( o) 


8z 2 


(17) 
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The resultant force and strain relations give 

e z (o) -4^ (o) - VK e (o) ] 


eg(°) = _i. 

hE 1_ 


9 


(o) _ lf ze ( ° ) 


hG 


Vjr z (°) 


( 18 ) 


The substitution of equations (l 8 ) and ( 16 ) into equation (17) leads to 


where 


v V°) = M 8 2 w (o > 

a 8 z 2 


7».2L + a -»L + 4 

4 2 2 k 

8 z 8 s 8 z 8 s 


s = a 9 


(19) 


The moment-equilibrium equation together with its relation to the changes 
in curvature (equations (B4, (B5), and (BIO) ) gives the second equation 

for w^°^: 


v V°) 


1 a g F (o) 

aD 8 Z 2 D 


( 20 ) 


Equations (19) and (20) are a pair of simultaneous equations for w^°^ 

and subjected to appropriate boundary conditions. The effect of 

curvature in the 9 direction is evident from the terms multiplied by 
l/a. Thus one obtains the flat-plate equations if (l/a) — ^ 0, that is, 
coupling of the two equations disappears. It is noted that these equa- 
tions are exact for bending of a circular cylindrical shell with small 
deflections. 
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• Equations (19) and. (20) can be combined into one equation for one 
complex function (cf. Reissner, reference 5). If equation (20) is 
multiplied by a constant k and added to equation (19) one obtains 
first 


then 


+ kw^ 


8z 2 


hE (o) kF(°) k P£ 

— w + — + — - 

a aD D 


V^[f ( o) + i^ED W (o ^] = i\/hEDw^ 0 i] + k (21) 

L J aD ^ z 2L d D 

if k = i\/hED where i = \f- 1. Equation (21) is a differential equation 

for one complex function + i\/hEDw^ °^J . Furthermore, if one desires, 

equations (19) and (20) may be reduced to one eighth-order differential 
for a function t such that 





v (o) 


a 

Eh 






(22) 


However, for practical calculation F^°) and w^ 0 ) may be represented by 
double Fourier series and equations (19) and (20) are used directly. 


Equations for Second Approximation 

The second approximation involves now the known functions of the 
first approximation on the right side of each equation. Let 

m (!) 

+ a(9,z) 


N 


> (1) ate.z) 


a 2 s© 2 


> 


N 


'©£ 


(!) = _i a 2 F (I) 


8 © dz 


(23) 
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where the superscript corresponds to this approximation. With equa- 
tions (23) force-equilibrium equations (B13) are satisfied if 


an 

ae 

an 

az 


2 || %z (0) + fN z (0) 


¥ - <°) + r - g H e <°> 


dz 


a z dz 


(2k) 


The relevant compatibility equation for this approximation is 


■\2 (1) 

a 

dz 2 


d 2 «. 


( 1 ) 


ae £ 


^ z e (1) 

a 0 dz 


= _1 
a 


a 2 v (1) 

az 2 


+ G ( °)(e,z) (25) 


where G^°^(0,z) is a function determined by the first approximation. 
With the aid of the resultant force-strain relations 


:u = jJh 

hEL 

:i) - ^ 
hEL 


.t 1 ) - VN U) 


■3 




( 1 ) - VN^ 1 ) 


L 3 


y U) _ _1_ N (1) 
7 z B hG a z6 


( 26 ) 


u 


and equations ( 23 ), compatibility equation ( 17 ) gives 


vM 1 ) = - M + n(o)( e ) 

a . Q ' 7 * 


a a z 2 


( 27 ) 


where H^°^(0,z) is a known function from the first approximation. The 
moment-equilibrium equation with equation ( 23 ) and the relations between 
resultant moment and change of curvature lead to 
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vV 1 ) = A ^ f(1) 
V ^ Sz 2 


+ K^(0,z) 


(28) 


where k(°)(0,z) is also a known function. The solution of equations (27) 
and (28)' would give the second approximation^ The effect of double 
curvature comes in at this stage. The method of solution follows the 
ones suggested previously. It is evident that the form of the equations 
for higher approximations will be similar in structure. Thus for the 
(n + l)th approximation one has 


V^p( n ) 

V 4 w( n ) 


_hE A 

8z 


+ h < 8 > z > 






Da a? 


(29) 


and K^ n-1 ^ contain all of the known (n - l)th func- 

w^ n-1 ^ and F( n-1 ^. The functions H^°^(0,z), g(°)(9,z), 

and K^ 0 Hs,z) are not given explicitly here but may be formulated 
readily from the expansions given in appendixes A and B. 


where 

tions 


Solution of Specific Problem 

The developed differential equations will now be applied to a 
particular problem. The problem at hand is that of bending of a curved 
plate with two meridian curves and two parallel circles as boundaries. 
This plate is loaded along the edges z = Constant in compression 
(i.e., the load is along the direction of the z-axis). The boundaries 
of the plate are given by z = ±Z and 9 = ±0 O . (See fig- 3* ) 

Boundary conditions .- Four boundary conditions are needed on each 
edge, and with four edges there are sixteen conditions for w and F 
together. 

The plate is supported at all edges and hinged at the edges so that 
bending moments are zero there. Thus, 

at z = tl : w = 0, M z = 0 

at 9 = ±6 0 : w = 0, Mg = 0 



(30a) 
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Furthermore, for applied load at edges z = Constant it is required 
that at z = tl 



d0 = Applied load 


(30b) 


(in most cases N z = Constant and the condition becomes equal to a 
constant. ) In addition, it is required that the shear resultant be zero 
at all edges, that is. 


at z = tl : N z g = 0 


at 9 = +9 0 : N ze = 0 


> 


(30c) 


There are two more conditions needed and these are supplied at the 
edges 0 = ±0 O . They may be formulated in either of two ways. If the 
plate has stiffeners attached at the edges which have infinite bending 
rigidity in the 0 direction then the boundary condition at 0 = i0 Q 
is 

u = 0 (304) 


and if the stiffeners have zero bending rigidity in the 0 direction 
the condition at 0 = ±0 q is 


N 0 = 0 


(30e) 


Of course, the actual condition would be between these two. It may be 
noted that condition (26c) may be modified to have zero tangential 
displacements at the edges (then N z g £ 0 at edges). With the help of 
expansions (12), ( 13 ), (1*0; and ( 15 ) the above boundary conditions 
imply the following conditions for the different approximations: 
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At z = tl : 


w = 0 gives =0 for n =' 0, 1, 2, . . . 

M z = 0 gives M^ 11 ) =0 for n = 0/ 1, 2, . . . 


At 6= +9 0 : 


w = 0 gives = 0 

% = 0 gives Mg^ = 0 for n = 0, 1, 2, . . . 

At z = tlx 

N z = Constant gives °) —Constant 
N z ^ =0 for n » 1, 2, . . . 


> ( 31 ) 


At e = +e Q : 

If N e = 0 then N 0 ^ = 0 for n = 0, 1, 2, . . . 

If u = 0 then =0 for n = 0, 1, 2, . . . 

At z = tl, 9 = ±9 0 i 

N z q = 0 gives N z g^ n ^ = 0 for n = 0, 1, 2, , . . 


Solution to first approximation. - For further application; consider- 
ation is given to stress-free edges at 0 = id 0 (i.e.; equation (30 e)). 

A possible solution to equations (23) and (27) is that 

v(°)= 0 


( 32 ) 
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and 


F (o) = Bs 2 (33) 

Equations (32) and (33) satisfy all boundary conditions given above 
and furthermore differential equations ( 19 ) and ( 20 ) are satisfied 
^with = 0^. The constant B is determined by the load at the 

edges z = tl. If the loads at these edges are: 


then 


N z = k = Constant 


> 


(34) 


B = k/2 


by equation (l 6 ). The displacement components u ^ 0 ^ and v^° ^ are given 
by 


» (o) - -1 - 


rlo> - Ef kz 


(35) 


Solution to second approximation .- The functions 
must first be determined from the first approximation. 


n, h(°), and K<°) 
Thus, one has 


= 0 

80 



> 


Sa = _i df k 

8 z a dz 


( 36 a) 


or 


n--££ 

a 


( 36 b) 
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H<o) (9 , z , = (2_^v)£ A + A £ 
a dz^ dz J a 


K (o) (e,z) 

dz D 


Differential equations ( 27 ) and (28) 1)6001116 now 


( 36 c) 


(36a.) 


vM 3 -) = -L a> 2 _F (1 ), .. x d 2 f 

Da 12^ D 17 


( 37 ) 


and 


^,(1) = + ( 2 - v )k d g f + d3f kz 


V F 


dz‘ 


dz3 


(38) 


Solution of equations (37) and (38) .- To satisfy the condition of 
simply supported edges assume a solution for w(^) in the form 


»(!).- 


cos 7 ®t z 


m=l, 3- • • n=l, 3. 


"mn 


■\2 1 


) cos (fg s ) (39) 


where b = a0 o . Insertion of equation (39) into equation ( 38 ) leads to 
V ^p( 1 ) =WL 5lI Vjftmi cos hn z cos v n B + 

m=l,3-*- n=l,3... 


(2 - V)k d 2 f d 3 f kz 
a dz 2 dz3 a 


(4o) 


where p m = mit/2 1 and V n = nrt/2b for abbreviation. The particular 
solution of equation (40) is 
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F p ( L ) = 


m=l,3. •• n=l, 3- 


F mn cos Sa z cos 


V n z + 


J(z) 


where 


F = _ 
mn 


h_E 

a 


7 2 b 2 

W* * v 2 )' 


W, 


mn 


> 


d^J(z) = (2 - v)k d 2 f d^f kz 
dz^ a dz^ dz^ a 


(4la) 


(bib) 


and 


= M- b 
'm m 


Si - v n^ 


A complementary solution that can satisfy the boundary conditions is 


(1) 


kf* j" 1 11 ■ 

= — cos M- z/A cosh u. s + C p. s sinh u s^ 

2a m= 13 * {_ m \ m m 111 111 m ) 

^ > J~cos v n s(B n cosh v Q z + D n v n z sinh v n zjj 

n=l,3... L J 


(42) 


To recapitulate, the boundary conditions which remain to be satisfied 
are, at z = +Z, 


a) s sV3.tt = 0 

z 3 b 2 


(*3a) 


m * 2 f (1) 

No 1 =-— = 0 

8z 8s 


(43b) 
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and, at s = lb. 




8 2 F^ _ 
Sz 2 

. aV 1 * 

5z 5s 


O' 


= 0 


(43c) 


(43d) 


Assertion of conditions (43a) and (43c) leads, respectively, to rela- 
tions between the constants of 


A* = 


V 


m 


cosh 7 


- 7 tanh 7 C 
'm 'm m 


m 


(44a) 


Bn = -S n tanh 5 n D n 


where f m is the coefficient of the Fourier cosine expansion of 
given by 


(44b) 

d 2 J(z) 

dz 2 


-p , ^ d 2 J(z) „ 

fyn = / COS (^X dZ 


m 


'-l dz £ 


(45) 


and 


Pm — 


nr 


a n = v n j 


The detail's of this calculation are given in appendix C. The remaining 
conditions, equations (43b) and (43d), give two infinite families of 
algebraic equations for the constants D n and C m . They are: 
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( Sn 2 + t*y 


sin sin cosh S n D n = 


- 7 m 3 b 2 
a 'm. 


n=l,3.\. ^ 

+ s “ 11 7 “) 


ol, sin £% ,2 

w mr , - — f m tanh 7 m - 


7 ^ + a 2 ) 
/ m n J 


2 "mn j 


for m = 1, 3, 5 


and 


Si 2 7m 2 Pm 


m=l,3... 

2*\ 2 
+ 7m) 

oo 

2 V7 ffi 2 sin 

z. 

m=l, 3. . . 



PmK 2 ' 

hEb 2 V 2- 


7m Pnfti 

a t~T 

m=l,3. 

1 • • 



^ sin Ojjj sin ^ cosh 7 m C m = 


m 


+ 7*r) 


- a n°n 


cosh 6- 


+ Sq sinh 5 nj 


n 


for n = 1, 3 ■ • • 


(b6) 


(Vr) 


The details of the above- calculation are also given in appendix C. At 
this point all boundary conditions are satisfied and it remains to 

satisfy differential equation (37)* If the value of is inserted 

into equation (37 ) > then this equation is satisfied if 

0 


w, 


mn 


,7m 2 + a n 2 ) + 


hb 1 ^ 

7 ^ 
'm 

- bVn - ^ 

a mn f m^ 

Da 2 

ro 

i 

iv) 

^ Dal 2a n 

P m 2 cosh 7 m 


2„ 2 


^7m a n 


4b* 


'm 


(“n 2 + 7 m 2 ) 


2 cosh 7 h sin On + 2 D n 


3 3s 2 

H m n 


^ (8n 2 ♦ S m 2 ) 


2 cosh 5 n sin P m 


for any m and n 


m 
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* 

where is the coefficient of the Fourier double cosine expansion 

of the function 


1 _ d 2 j(z) d 2 f k 
D dz 2 “ dz D 


and is given by 


Kmn 



l fd 2 J(z) 

DL dz 



cos ji^z cos v n s ds dz 


( 49 ) 


The complete solution of the problem lies in, the solution, of the infinite 
families of algebraic equations for C m and D n and then v mr) from 

equations ( 44 ). An approximate numerical solution for C m and D n may 

be obtained as follows. It is assumed that coefficients with suffixes 
greater than some fixed number can be neglected. Then it is verified 
that increasing this number does not affect the coefficients with small 
suffixes. 

Illustrative example .- As an example the following meridian curve 
is considered: 


r = a 



(50) 


It is easily verified that the slope dr/dz is zero at z = il . Thus 
the load N z = k at z = +Z is normal to the edges. In addition; con- 
sider a plate such that Z = b. The sign of the constant k determines 
compression or tension; that iB, k < 0 indicates compression and 
k > 0; tension. The procedure in the numerical calculations is outlined 
in appendix D. The results of the calculations are summarized in table I 
and figure 4 . It may be noted that all the calculations can be carried 
through with the sign of k arbitrary. Thus the sign of k which 
determines either compression or tension can be assigned at the last 
stage of the calculations. 


DISCUSSION OF RESULTS 


The numerical results for the example are given in figure 4 . It 
is seen from this figure that maximum value of the deflection w 
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(normal to undeformed plate surface) is not at the center of the plate. 
This deflection appears to have a maximum value near the loaded edge. 

This is due probably to the choice of the form of the meridian curve 

which gives a negative Gaussian curvature for z > l and a positive 

Gaussian curvature for z <-— l . The deflections are largest at the 

{3 

part of the plate with negative Gaussian curvature. The moment distri- 
bution curve shows the same trend. It is noted that the actual values 
of the deflections and moments are multiplied by the parameter p. The 
order of magnitude of the parameter p is about 0.01 to 0.02. Thus 
the order of magnitude of the normal deflection w is about 0.04 to 0.08 
of the thickness of the shell. This indicates that the first two approxi- 
mations of most problems (i.e., p° and p^-) will be sufficient to give 
a complete solution. 

The convergence of the approximate solution of the infinite families 
of algebraic equations is quite rapid. These results are given in 
table I. A total of six coefficients was used; the dominant coefficients 
were the first two, and D]_. 


CONCLUDING REMARKS 


A method has been developed for -the analysis of the deformation of 
doubly curved thin plates under edge and surface loads. Only small 
deflections (small compared with the thickness of the plate) are con- 
sidered here. This method is particularly suited to the analysis of a 
plate with a large radius of curvature in one direction. This is clearly 
the situation existing in airplane fuselages. 

For the problem of a doubly curved plate in edgewise compression 
two infinite families of algebraic equations were obtained in order to 
satisfy the boundary conditions that would exist in most airplane 
coverings. Results were obtained by replacing these infinite sets by 
finite ones (neglecting all coefficients beyond a certain suffix) . Con- 
vergence of this approximate solution of the algebraic equations was 
quite rapid. 

Since the order of magnitude of the parameter p is small (in the 
hundredths) the first two approximations are usually sufficient to give 
the complete solution to similar problems. Moreover, the convergence 
of most expansion methods about a parameter is difficult to prove mathe- 
matically. However, in most practical applications the parameter p will 
be small and only a few terms in the expansion will suffice for the solu- 
tion. Of course, the solution so obtained will indicate the number of 
terms required in the expansion. 
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It should he noted again that the equations are valid only for 
deflections that are small compared vith the thickness of the plate or 
shell considered. If, in the application of the method, large magnitudes 
.of the deflection are obtained one must resort to a nonlinear theory. 

A specific example is presented which . can be immediately applied 
to the bending of thin plates with compound curvature for small deflec- 
tions when loaded in edgewise compression or' tension. The basic equa- 
tions given are equally valid for a curved plate loaded by any other 
edge loads or by surface loads, but the form of the specific example 
as presented will change slightly. Moreover, any shell of revolution 
with a meridian of small curvature may be analyzed by this method. 


The Pennsylvania State College 

State College, Pa., August 28, 1951 
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APPENDIX A 

EXPANSIONS OF TEEMS IN POWERS OF PARAMETER p 


To facilitate computation the following expansions are particularly 
useful. All terms are expanded up to and including second powers in the 
parameter p. With 


r = a + pf(z) 

The following expansions can "be immediately written: 


1 _ _ _1 
R© “a 


f 2 


f 2 


1 - u. — + - — - «=• f - 

^ a M a 2 2 z 


(Al) 


— = uf 

Rg z 


a = 1 + p 


2 V 


§7 = ^ af zz + ^ 2ff zz 


1 

r 




i = 1 - p 2 


-z _ 


_ „2 f 


-Z = p -*• - p- -% r z 
ra a a 2 z 


Differentiation with respect to z is indicated by the subscript z. 
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APPENDIX B 

EXPANSIONS OF EQUATIONS 
Expansions of Equations for Strains 


The strains have the following expansions up to and including the 
second powers of ,.p: 


(o) 

_ i a«‘»> 

(0) 
w v ' 





a 80 

a 




(1) 

_ 1 8u^ 


_ x 

8u^°) 

+ £z v (o) 


a 

a 

a2 

80 

a 

(2) 

_ l 8u^ 

(2) 
w' ‘ 

_ t_ 

8u ^ 1 ^ 

+ £2 v (1) 


a 80 

a 

a 2 

80 

a 


f 2 

ff Z 

y( °) 

.11 

_ I f 2\w 


a^ 80 

' a2 


U 2 

2 f z J • 


\ (Bl) 


(o) _ 8v^°) 

2 ‘ 8 Z 

(!) _ $v^ + f w (o) 

5 8z zz 

(2) _ dv( 2 ) , f Jl) 1 _ 2 8v(°) 

2 - IT" + fzzV - '-5 fz "87” 


> 


(B2) 
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( 0 ) 

du (o > 

1 

dv^ 





?ze 

dz 

+ — “ 
a 

d0 





y (1) 

du^ 

+ 1 

^v(l) 

£z 

u (°) 

^ X 

dv ( °) 

'z 6 

dz 

a 

ae 

a 


a 2 

de 

' i 

/ 

7 S 2 } 

_ du^ 

+ I 

av< 2 > 

£z 

u^) 

_ x 

av (1) + 

/ z 9 

dz 

a 

de 

a 


a 2 

ae 


ffz U 

( 0 ) 

f z 2 


, f 2 

dv(°) 



2 

dz 

a 3 

6e 

■ 


( 0 ) = - -i. 


1 aV o) 


£ 3a2 


«e (1) 


1 

‘"2 

a 


a 2 ™' 1 ’ ^ 2f a 2 >> 


f * a>> 


ae £ 


a^ d0 2 


a 


dz 


K (2) 1 d 2 ;/ 2 ) , 2f d 2 ^ 1 ) fz ch^ 1 ) 

Kq ' * = m. — ^ ^ T ~ ?S — T ' 


a 2 d0 2 


i3 ~W 


dz 


(Bh) 


2 f\d 2 v^ 


J- 2 ±_ + i, 

a 2\ a 2 a y d©2 


J 


K (o) _ _ 8gy..?l 

Z “ dz 2 
K (l) = - ^ - 

2 8z2 


K (2) _ S 2 ^ (2) _ + 1 f 2 

KZ dz 2 2 Z 


a^o) 

Sz2 


d>> 


+ fzfzz 


£z£ 5^°) 
2 dz 2 


> (B5) 
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T (°) = 2 8 2 y(°) 

" a 30 3z 

t ( l) _ _ 2 3 2 w^) + £z 8y(°) + _f_ + £z 8y(°) + _f_ 8 2 w~( 0 ) 

a 3z 50 a 2 3© a 2 80 3z a 2 80 a 2 8z 50 


t < 2 ) = 


.2 8 2 v^ 2 ^ 

a 5z 80 


i 2 5z 80 


f z , f 8 2 w^ , f z 8v^ 1 ^ 


1 2 80 


a 2 3z 50 




80 


- ^ (B6) 


2ff„ 


8w(°) 

80 


• z 
a 


8^°) 

50 8z 


18 c 
a^ 8z _ 


8w^°^ j.2 

80 


f. 5^°) 
3.3 8z 50 


f z g 8^°) 
2a 8z 80 


Compatibility Equations 

The following compatibility equations can be constructed from the 
values of the strains given above: 


a% (o) , 1 &V o) 1 a 2 7 z8 (o) 1 a^o) 

8z 2 + a 2 50 2 ‘ a 8z 50 " a 8z 2 

aVu ^ T &,(u x a\ fl d) , ft 2 w( i) f ,- g2>) 

8z 2 a 2 80 2 a 8z 50 --- a 8z 2 a^ 80 2 

f zz (o) . £z 8y(°) , £z 8w^°^ , X d 2 ^ 0 ) 
a 2 a 2 8z a 2 8z a 2 5z 2 

£z 8^°) f 8^u^ 0 ) + 2 ^zz 8v^°^ + £z 8 2 y( 0 ) + 

a 2 50 8z a 2 30 3z 2 : a 3z • a 3z 2 


(B7) 


+ 


f zzz v (o) + X f 83 y(°) + £z 8 2 v(°) 
a a 3 30 2 3 z a 3 30 2 


(B8) 
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The compatibility equations for the next and succeeding approximations 
are constructed similar to equations (B7) and (B3) obtaining in each case 
, n 2 (n-1) 

terms as - - — ^ and so forth where n denotes the nth approximation. 


Expansion of Equations for Moment Equilibrium 

The last three equilibrium equations (equations (4c), (5a), and 
(5b)) can be combined into the single equation 


a + 1 d 2 (iM z e) , 1 ^z£ + 5_ 

r 50 2 r 5z 50 r 5z 50 5z 


isK) 


5%i 


z 0 


50 5z 




(B9) 


The expansions in powers of \i of equations (b 6) lead to the following 
equations for each power of \i: 

For |u°: 


1 a W o) . a W o) , a V o) a W o) 


a 50 2 50 5z 


0^02 ( 0 ) 

a — + __Z| + + p.a = 0 

5z 2 50 5z 6 i 


(BIO) 


For m- : 


1 a V 1} + + a ^z (I) + ^0z (1) + N (1) + f 

a 50 2 5z Sz 2 50 5z 6 p £ f " 




5m. 


f _ 9 2 - ^zg 

aF 50 2 a Z 50 


(o) 


5m 


(o) 


+ - fz — ^ + 2f z + f zz ^ o) + 


5z 


5% (°) 


f h Z tl - fz ~Y Z f z&e KUJ - af zz Nz v 


- af zz N z (°) = 0 


(Bll) 
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For 


P 


2 . 


aV 2> + a \ z (2) + a 3 2 m z (2> + (2) 

c) 0 2 50 5z 5z 2 5S 5z 


z f zz 2a 

+ + p — + 


■e 2 


( 2 ) _ ^ w ( 1 ) _ f_ ^ 2 M e (I) + f 52 m z (1) + 2 £z ^ez (1) 


V~' - af zz N z 


a 2 as 2 


bz 


S0 


5m 7 ^ ( 1 ) 

2 f„ -4 + £„JL K ' 


z 5: 


L zz“z 


3m„ (1) a) 

f z -fj- - *£iV ' - ff 


n (») ; 


zz^’z 


f z 2 a 2 ^ z (o) _ , ^ 0z ( o) 


+ ± 1 z _ a. f 

2 a 2 / a 5 © 2 -2 z 5z 2 


- af 7 f 77 — 

Z ZZ 5z 2 


= 0 (B12) 


Expansion of Equations for Force Equilibrium 


The remaining force-equilibrium equations may be expanded in a 

s imi lar way and then if each factor of p°, p-*-, and p 2 p s equated to 

zero the following equations result: 


For p°: 


a 50 

5u z (°) 

5z 


Sn z6 (0) 

+ — — = 0 

5z 

1 an z s (o) 
a 50 


> 


(B 13 ) 


For p 




H 6 
a 50 

5N J 1 ) 


z 0 
ciz 


+ 2 


%K z e ( ° ) * f 


, 1 aa z9 

a 


( 1 ) 


50 


+ -a N z 
a * 


( 0 ) + £ !2& 


an ze (o) 

5z 

(o) 


'■’N 


= 0 


> 


(Blit) 


5z 


N0 


(0) = 
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For |i : 


a% (2) + aN z 0 (2) 


ae 


az 


r z a) f aN z e^ f z 2 a%(°) 

+ 2 — N ze ^ ; + - — 52 + — = 0 


a Sz 


2a 


ans (g) + 1 aN ?9 ( 2 > + fz h jd + f W 1 ) _ £z % <i) + 


az a 5$ 


a 5z 


1 f z 2 a « z e (o) 

2 a ae 


= 0 


(B15) 
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APPENDIX C 


EVALUATION OF CONSTANTS OF INTEGRATION 


The expression for the stress function F chosen to satisfy the 
differential equation for F and the boundary conditions is 

F (l) = £ Z! F mn cos ^m z cos v n s + J ( z > + |f ® 2 + 
m=l,3 n=l,3 

00 °° 

Z cos p m zF m (s) + cos V n BF n (z) (Cl) 

m— 1 y 3 » • • 11=1 j 3 * • • 


where 

F m (s) = Aja cosh ^s + C^ m s sinh 


F (z) = B cosh V z + D V z sinh V z 
n' n n n n n 


and ^ = mrt/22, V n = rift/Zb, and s = a0 for abbreviation. Equa- 
tion (Cl) satisfies differential equation (21) if 


d^J(z) , 2 vk d 2 f t kz d 3 f 

dz^ a dz^ dz 3 


(C2) 


The boundary conditions which must be satisfied by F (equa- 
tion Cl) are, at z - tl: 


„ (1) . aVil „ kf = 0 


8s 


(C3a) 


N (1 ) _ 8 2 f (1) 
W Z0 - “ZTT7 - 


= 0 


(C3b) 


dz 8s 
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and, at s = ±b: 


% 


(1) _ s 2 ? (1) _ ? 

3z 2 


N 


zQ 


(1) = ■ ^ 2 f (1) = 0 


3z 3s 


Conditions (C3a) and (C3c) above lead to relations 
of constants and C m and B n and D n . Assertion 

tions leads to 


CO 



m — 1 y 3 • • • 


hn 2 F m (*) 


cos ti m z 


d 



and 

CO 

) cos v n B = 0 

n=l, 3- • • 


Equation (C5) immediately gives 


B n = - D n 8 n S n 


where 5 n = V n Z. 

Equation (C4) is a Fourier cosine expansion of the 
d 2 J(z)/dz 2 and its coefficients lead to 


Am = 


‘1 



cosh b 


p m b sinh i^b 
cosh p m b 


Cm 


(C3c) 


(C3d) 


between the sets 
of these condi- 


(c4) 


(C5) 


(c6) 


function 


(C7J 



36 


MCA TN 2782 


where 


m 



^ ^X Z -^- cos i-^z dz 
dz^ 


(C8) 


To satisfy conditions (C3b) and (C3d) the following relations from 
equations (Cl) and (C3) are obtained: 


Mm sin ^ l F m '(s) = ZI 
m=l , 3 n-1, 3» • • 



F mn^m v m sin ) " 


v n V( J > 


sin v n s 


(C9) 


and 


n=l, 3. • • 


V n sin V n b F n '(z) = 21 

m=l,3. ■ 


2_ F mn v n^m sln v n b ) " 

l n=l , 3 • • • 


hn F m' 


sin p ffi z 


(CIO) 


where the prime indicates differentiation .with respect to either z 
or s. Equations (C9) and (CIO) are Fourier sine expansions of func- 
tions F m '(s) and F n '(z), respectively. The coefficients of these 

series with the help of equations (C6) and (C7) reduce to the following 
sets of infinite- algebraic equations: 
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m=l, 3 . • • 


4 (3 m a m 2 >'m 2 7 . , „ 

— — cosh 7m sin sin C m 

\, a n + ?m / 


2 V 7 m 2 sin “n sin P m 


m-i^ 3. P m (“ n 2 + 7 m 2 ) 


+ Z_ F mnPm a n sin P m “ 

m=l, 3. . . 


° 4 i 


5 n 

cosh 8 t 


+ 6 n sinh 5 ^ |D 


n n 


for n = 1 , 3 , 5 - 


(Cll) 


and 


CO 

Z 


2 q 2 
’m 


n=l, 3 . . . 


SiV f> 

l( 8 n 2 + Pm 2 ) 


2 cosh sln Pm 8in “n D n 


00 



n=l, 3 . . . 


WmSi 8in a n - 


tanh 

.Pm 


^m^m! 


r m 


cosh 7. 


+ sinh tJC 


m/~m 


m 


for m = 1, 3,5. . . 


(C 12 ) 


■wh ere 3 m = P- m l, = v n b, 5 n = v n Z, and 7 m = p m b. There exist also 
two relations given in the text whereby may be eliminated so that 

the two sets of equations contain only the unknowns C m and D n . These 

two infinite families of algebraic equations may then be solved 
approximately. 
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APPENDIX D 

PROCEDURE IN NUMERICAL CALCULATIONS 


The procedure in the numerical calculation follows here. First, 
equations (42), (43), and (44) are made dimensionless by using the 
quantities 


- CL 
r>_ = m 

ki d 


W = 

w mn h 


- £ 

kZ‘ 


n - n 

Dn _ __ 


5 - X™* 0 ^ 

nnn - — — ' 


f = !e 

m ki 


(Dl) 


Then equations (42), (43) , and (44) become 


WW = Kmn + S mn + C m T mn + D n V mn for ^ ( D2 ) 


t. ^-W)- £ + '£ 4b* 

m=l, 3- • • ' «mn / m =l,3... m=l,3... \^mn 




^mn , 


DnfZn + b 

\ m=l,3... 


iffln 


T mn 

^mn. 


for .n = 1,3,5- 


(D3) 


£ II Jmn(^^ ; 

n=l,3« • • ' %n / n=l,3... \ Kffln 


Dm + Cja(M m + 


— T \ 

p mn 

n=l, 3* 


(D4) 


■'mn 
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These functions are tabulated, at the end of this section. The physical 
parameters which appear can be combined into three groups: 

kb^ h 2 Eb 2 
haD 

As a numerical example, assume that 


> (D6) 


These values correspond, for example, to a fuselage panel with these 
dimensions: 


£ = 600 
h 


kb" 


= 10 


hEb: 

Da 2 


(D5) 


h = 0.03 in. 
a = 60 in. > 


b = l = 18 in. 


(D7) 


The procedure in the numerical example 
(l) The particular solution J(z) 


is as follows: 
is given by 


J(z) = 4- — - ( 2 
30 


(2) The Fourier coefficients 



1 _1_ + ^( 2V ~ jj 1 + 1 - V 2k 

(2 - v) P m (2 - V) p m 3 2 - V p m 5 


f m = -2k(2 -v)l sin p m 
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(3) The Fourier coefficients are: 


K mn DbZ 


r b r 1 

d 2 J(z) d 2 f v 

1 -b J -l 

dz 2 dz 2 


cos p m z cos V n s ds dz 


2k sin sin 

Da cx„p_ 
n r m 





k + 



(4) The functions as defined "by equations (D2), (D3), and (D4) are 
tabulated and summed where required. Values of these functions are 
given in table II. 

(5) The two families of algebraic equations are then set up for a 
finite number of unknown coefficients D n and C m neglecting all coef- 
ficients beyond a certain suffix. (in the numerical example a total of 
six coefficients was used.) These equations are then solved by Crout's 
method. The results are given in table I. 

The functions used in the numerical example are defined as follows: 


^mn 



+ trf + ^ 
^ ' Da 2 

kb V haD 



Kmn = K^aD/k 


mn 


2f ffl a n sin a n 

2 2 
a n + 7m 
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T = 
mn 


^m 2ct n^m C03h Bln a n 
( a n 2 + 7 m 2 ) 2 


V. 


mn 


_ 


(5n 2 - Pm 2V 


75- cosh 5 n sin 0. 


m 


w mn = : g sin P m cosh ? m sin 


(“n 2 + 7 m 2 ) 


^mn _ “ 


2 Vm sin gn sin 

2 2 
“n + 7 m ^ 


h 2 Eh 2 7m 3a n 


Kft (r m 2 + vf 


si^ Pm 


Z„ = - 


n \cosh 5 n n 


+ 5„ sinh 5 


n , 


4s 2 & 3 

I mn = 7 — : m n _ Q sin “n sin P m cosh 5 n 


(5n 2 + Pm 2 )‘ 


T _ hV 7 i 3 “c eln <% 

J mn o o 

kl^a / 2 2 \ 

Vm + “n ) 
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TABLE I 

RESULTS OF NUMERICAL EXAMPLE 


(a) Deflection coefficients and values of 
moments at center line 


(Upper sign 


for compression and lover sign for tension 
vhen tvo signs are present] 


z /l 

MQ^b^Dli 

v^Vb 

s/b 

M z ^ 1 ^b 2 /Dh 

0 

*2. 4191 

±3.00997 

0 

±25.432 

.2 

±10.374 

± 2.6122 

.2 

*17.892 

• 5 

±30.767 

±3.4933 

• 5 

±3.6157 

.8 

±71.944 

±3.85909 

.8 

±114.240 

1.0 

0 

0 

1.0 

0 


(b) Values of unknovns for different 
approximations of algebraic equations 



Number of unknovns 


2 

4 

6 

Cl 

-60.384 

- 52.232 

-51.464 

d 3 


-.00265 

-.005040 

*5 



. 000100 


- 1.6976 

-1.70914 

- 1.7547 

°3 


. 049589 

. 049369 
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Figure 4. - Deflection normal to shell surface and moment distribution 

at center lines . 
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